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SIGNAL ANALYSIS AND ERROR ANALYSIS STUDIES 
FOR A GEOPOTENTIAL RESEARCH MISSION (GRM) 

W. D. Kahn 
T. L. Felsentreger 

ABSTRACT 

The Geopotential Research Mission (GRM), formerly known as GRAY SAT/M AGS AT, has been 
proposed for the mapping of the Earth’s gravity field, It is a candidate for a NASA new start in the 
near future. GRM is composed of a pair of surface-force compensated satellites in identical polar 
orbits, spaced approximately 3° apart at an altitude of 160km with a system precision of I 

The objectives of this paper are to discuss the signal characteristics and the geopotential param- 
eter recovery capability of the SST Doppler sensor flown on GRM. 

Simulation studies of the velocity jvrofiles resulting from the perturbation produced by a 1° X 
1°, 1 mgal anomaly as sensed by tv/o GRM spacecraft orbiting altitudes of 160km and 200km re- 
spectively are described. It has been found that the amplitude of the gravity signal drops off by a 
factor of 1.5 when going from an altitude of 160km to 200km. By extrapolation the signal ampli- 
tude is further decreased by a factor of 3 when the orbital altitude is increased to 250 km. Thus the 
amplitude of the measurement drops off as the altitude is increased to the point where it is insignifi- 
cant at the 1 mgal level for altitudes above 200km. 

Spectral analysis results show that for a GRM mission altitude of 160km and a system precision 
of 1 micrometer/sec., gravity field information can be sensed up to 230 cycles per orbital revolution - 
beyond that frequency the gravity signal is characterized by white noise. It follows that at the GRM 
mission altitude of 160km and a satellite-to-satellite Doppler system precision of ±1 micrometer 
per second, 1° X 1° gravity and geoid anomalies can be determined to an accuracy of 3.4 mgals and 
8.6cm respectively. 

These results are compatible with the scientific requirements for GRM which are; Accuracy - 
2.5 mgals; Horizontal Resolution - 100km; Accuracy of Oceanic Geoid Undulation Difference - 
lOcm (over wavelengths of 100 to 3000 kilometers). 
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I. INTRODUCTION 


Studies have been conducted at the Goddard Space Flight Center to investigate the accuracy 
with which the Earth’s gravity field can be determined from a gravity field mapping mission such as 
the Geopotential Research Mission (GRM), fomierly designated as GRAVSAT/MAGSAT. The 
studies performed at Goddard are threefold: (1) gravity field signal analysis, (2) rapid error analysis 
which can be used for mean gravity anomaly and mean undulation accuracy estimation for almost 
continuous data, and (3) full-scale simulation studies, Jhis paper reports on the results of (1) «& 
( 2 ). 

II. ANALYSIS 


A. Signal Analysis 

The perturbation by the gravity field induced on the intersatellite Doppler signal is given by 
the following relationship: 


rh ^ 

Hrh = r 


1 

r d\p 


dt 
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dt 


(1.0) 


where 

T = disturbing potential 
r = a + h 

a = Earth’s mean equatorial radius 
h = satellite orbit altitude 
\jj = true anomaly 

5t = time separation interval of two satellites in identical orbits. 
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The form of tlie disturbing potential most useful for formulating ( 1 ,0) is given by the formula 
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Ag = mean gravity anomaly 
do = element of area. 

Using (1 .0) and (1,1) for an orbit configuration as depicted in Figure 1, velocity profiles resulting 
from the perturbation produced by a I mgal gravity anomaly as sensed by the 2 satellites a"e shown 
in Figure 2, The satellites are spaced 330 km (3°) apart ~ results are shown for 160 km and 200 km 
satellite altitudes. The amplitude of the gravity signal drops off by a factor of 1 .5 when going from 
an altitude of 160 km to 200 km. By extrapolation, the amplitude is further decreased by about a 
factor of 3 when the orbital altitude is increased to 250 km. Thus, the amplitude of the measure- 
ment drops off as the altitude is increased to the point where it is insignificant at the 1 mgal level for 
altitudes above 200 km. By way of contrast. Figure 2A shows the SST signal generated by GRM for 
one full orbit, under the influence of a full (36, 36) spherical harmonic gravity field (GEMIOB). 

The two satellites are at 160 km altitude, separated by 3 degrees. 

The ability to resolve the signal from a (1° x 1°) 1 mgal gravity anomaly for 2 spacecraft orbit- 
ing at a 160 km altitude, whose orbits are displaced from the anomaly by 1“ and 5°, respectively, is 
shown in Figure 3, From a 1° orbit displacement, better than 60 percent of the gravity signal is 
sensed. However, when the orbit is displaced by 5° from the gravity anomaly, only 1 0 percent of 
the gravity signal- is sensed. This would imply that the gravity anomaly signal from an adjacent 
block drops off drastically with increased distance. Nevertheless, it is significant that remote gravity 
anomalies do contribute to the overall gravity signal by as much as the 10 percent indicated. 
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B, Error Analysis 

A rapid error analysis procedure developed by C, Jekeli and R. H. Rapp (Ref, 1) has been used 
for the estimation of mean gravity anomalies uSul mean undulation accuracies that can be obtained 
from SST missions such as the future Geopotential Research Mission (GRM), The basic method 
used is the analysis of the gravity field mapping in the frequency domain with results obtained in 
tenns of mean gravity anomalies and mean geoid undulations. The measurements to be considered 
arc velocity differences, with the errors in the measurements assumed to be uncorrelated. 

In order to estimate the expected "signal” for SST Doppler we shall evaluate its power spec- 
trum (Refs. 2, 3) derived from the Earth’s anomalous gravity potential. Tliat is, 

(V+T) = — [1+ J I (a/r)® P 5 j„(sin^)(Cg„^cosmA + Sg„, sinmX)) (2.0) 

I* *0 

where 

r is the geocentric distance 

if is the subsatellite geocentric latitude 

X is the siibsatellite east longitude 

a is the mean radius of the Earth 

Pjgjj^(sin ip) are the fully nonnalized associated Legendre polynionials of degree and 

order m (£ is also the orbital frequency of a satellite in cycles per revolution) 
and 

Cgni> S arc the fully normalized coefficients of a tmneated harmonic series expansion 

representing the geopotential. 

The power spectrum of the velocity can be determined by utilizing the Law of Conservation of 
energy, 

>/2V^ - U = E (3.0) 

where U is the total gravitational potential, v is the velocity, and E is the total (constant) energy. 
Decomposing U into the reference potential V and tlie disturbing potential T of the Earth, we have 
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i/jv2 ^ V - T « E 


The variation of (3.2) then reads 


v6v « 5T 


(3.1) 


(4.1) 


where - GM/r. 

From Equation (2.0), together with (4.1) we obtain the variation of the r udUte velocity mag- 
nitude as: 

5v » V ^0 ^J2m 

The velocity power spectrum fjf 6v is now obtamed by squaring Equation (4.2) and integrating 
over the unit sphere: 

V/(5v) = v2(a/r)2^ Pfim (sinv?) (Cg^ + (4.3) 

Using Kaula’s rule (Ref. 4), the summation containing the C, S terms can be estimated as follows: 

,i (^L + SL) = (4.4) 

Then 


(5v) = v2 {mlxf^Cl9. + 1) (10-^12^)^ 
> 2 


(4.5) 


In order to derive the velocity variations for the low-low SST configuration, Equation (4.1) is 
to be modified as follows: 

5T = (T 2 ~ T^) = v5vh (4.6) 

where Tj^ and Tj are the disturbance potentials at satellite positions 1 and 2 in the orbit with 
an angular separation, i//j 2 - 
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Tlie corresponding vevtonce in the horizontal component of velocity is now to be derived as 
follows: 

(6T) = <7j 2(T,) + aj^CTj) - 208 (T,Tj) (4,7) 

By assuming that = “S^CTa) = fg^CT) and ujCTiTj) = Pj;(cos i^'u) oj®(T) (Ref2) 

one obtains 


oj 2(ST) = 2aj,*(T) [1 -Pj(cos(tjPl (4.8) 

and from Equation (2,0) together with (4.6) and (4,8), 

ag2(6 a^^) - 2v^ (a/r)^^ 1 1 - P£(cos (22 + 1) (10“^/£2)2 (4 9) 

Equation 4.9 represents the horizontal velocity “signal”, variance. 

The intersatellite velocity signal spectra is given in Figure 4. It can be seen that for the planned 
GRM altitude of 160 km and the SST system precision of ±l Mm/s, gravity signals at the orbital 
frequency of 1 80 cycles/revolution (2 = 1 80) can be sensed by GRM. Also shown in tliis figure is 
the level of signal attenuation with increasing orbital altitude. 

Tlie frequency spectrum only provides information on the specific frequency (spectral line) at 
which the SST velocity signal will sense the gravity field. It does not sense the maximum frequency 
at which a SST system with a specified signal precision and a specified orbit altitude can sense 
gravity field information. Determining the total gravity signal sensitivity for a SST system having a 
given level of precision necessitates integrating the velocity cpectrum. That is. 


o^2,co(5v„) = 



(5.0) 


The lower integration limit can be interpreted as that orbital frequency at which the gravity 
signal-to-noise ratio is unity for a given SST system precision. In Figure 5, it can be seen that for 
the system precision of 1 /im/sec and orbital altitude of 160 km, the gravity signal cut-off frequency 
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occurs at the orbital frequency of about 230 cycles per revolution (8 = 230). All gravity iiiforma- 
tion above that frequency (i.e. limited by the SST sensor precision) would be characterized by 
white noise. 'IHius, a higher “cut-off’ frequency at a given orbit altitude necessitates a SST sensor 
of a higher precision. 

To assess the accumey with which geophysical parameters (such as gravity or geoid anomalies) 
can be determined when using the previously selected GRM sensor precision of ±1 ixmls, the models 
developed by Jckeli and Rapp (Ref. 1) are used. 

The total error budget associated with these parameters is the combination of the errors of 
commission and that a omission (generally the larger terms). The models for these error sources 
read as follows: 

(a) The error of commission in terms of mean gravity anomalies is given by: 

f S, /3p2(S-l)2 (22+1) /r\2«+n''^ 
Oo(6S) = (7/a)« (AA/4,)« c(6v„) [ W J 

(b) The eiTor of omission in terms of mean gravity anomalies is given by; 


, r 0 (2~ 1) (22+ I)! 

a„(6g) =. TX 10-s I Ps,^ i 


where 


^ cot to 

2 L 2(2+1) J 

T = (GM/a^) 


2j = Value obtained from Figure 5 for which the signal-to-noise ratio is unity for a specified 
SST system precision ct (6v„). 

AA = Area of block in units of radians squared. 

\pQ = The radius of a circular cap that has about the same area AA as a block whose sides are 
S° (in degrees). 
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Angular separation of the two GRM spacecraft (in degrees), 
a =s Mean equatorial radius of the Earth. 

S® = Side of anomaly block (in degrees), 
and 

(c) The total enor for mean gravity anomalies is as follows: 


aj.(5g) = [a^'^(6s) + (^o^(5g)]''^ 


(6.5) 


In terms of geoid anomaly errors, the error models arc as follows, 

(a) The error of commission in terms of mean geoid anomalies is given by: 


ffj,(N) = (a/ 7 ) 
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rSi (Jj»(2S+l) 

(AAMt)** o(6v„)J „S - — — -—I — 

"^12=2 2[I-Po(cos i//io)]\ a/ J 


-P£(cos 

(b) The error of omission in terms of mean geoid anomalies is given by: 


OO 

Oo(N) = ax 10-5 [ S |3g2 (2£+l /£-»)] '^ 
1*^*1 


(c) The total error for mean geoid anomalies is: 


(7.0) 


(7.1) 


= [ac^(N) + 0^2 (N)] 


■A 


(7.2) 


Figures 6 and 7 show the total expected errors in terms of the gravity and geoid anomaly 
parameterization as obtained from Equations (6.5) and (7,2), The errors shown are for the GRM 
mission SST sensor precision a(5vjj) = ±1 fim/s, satellite separation distance ~ 3° as a func- 
tion of orbital altitude and block size. 

It is shown that by using a pair of spacecraft in a 1 60 km polar orbit separated by 3°, one can 
map the geopotential in terms of 1“ x 1° blocks to about 3H milligals or P cm respectively. 


CONCLUSIONS 

As is evident from Figures 6 and 7, the planned GRM mission requirements can be met using 
the SST low-low approach and a 1 jiim/s Doppler system, that is, global gravity and geoid anomalies 
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Figure 6, Gravity Anomaly Error (mgals). Figure 7* Geoid Undulation Error (cm). 



can be recovorotl witli a resolution of 100 km ami a precision of 3 msals and 9 cm 
respectively. 
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